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1 Introduction 

Let if be a number field containing the field of Gaussian rationals C K . fn 
these pages we are mainly interested in the computation of if— rational points of 
zero-dimensional algebraic varieties given by systems of multivariate polynomial 
equations. Namely, let /i, . . . , /s G . . . , X„] be a sequence of multivariate 

polynomials with integer coefficients. Let V{fi, . . . , fs) C C" be the complex 
algebraic variety of their common zeros, i.e. 

V{h, ...,/,):= {a; e C" : ^0,l<t<s}. 

For sake of simplicity, let us assume that V{fi, . . . , /s) is a finite set (i.e. a zero- 
dimensional algebraic variety). The set of if— rational points in V{fi, ■ ■ . , fs) is 
the set of common zeros of the system /i, . . . , /s whose coordinates lie in if", 
namely 

VK{fi,...Js):={x€K^ : h{x) = ■ ■ ■ = fs{x) ^ 0}. 

The goal of these pages will be to discuss several aspects of procedures perform- 
ing the following task : Assume that the field K is fixed. Given the sequence 
/i, . . . , /s, compute all if —rational points in VkHi, ■ ■ ■ , fs) (or eventually all 
if —rational points in Vxifi, ■ ■ ■ , fs) of bounded height). 

Note that the assumption on the field if is not very restrictive : For every 
zero C. G ^(/i, ■ ■ • , fs), there exists a minimal number field if (C) containing all 
the coordinates of C- The degree of the field extension if (C) over Q can also 
be denoted by deg(Q. In the sequel, the degree [if : Q] may be replaced by 
deg(C), and the results will equally hold. 

For our study, we consider a precomputation task which prepares the input 
F := (/i, . . . , /s), before we study the desired if —rational points. Procedures 
performing this precomputation task are usually called multivariate polynomial 
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system solvers applied to the input -F. The output of such polynomial system 
solvers is called the solution of the system F. Observe that all usual notions 
of solution of F will yield a description of the variety V{fi, . . . , fs) (cf. also 
[|CGH+99[). 



Here, we consider two (conceptually different) notions which define what a so- 
lution of the system F should be : coming from different fields, the notions are 
related to a symbolic/geometric and a numerical analysis/diophantine approxi- 
mation context : Kronecker's geometric solution and Newton's approximate zero 
solution. 

Thus, our study includes a comparative study of both approaches with regards 
to the basic problem described above. It must be said that our study is not 
intended to be either complete or definitive. It just tries to point out some 
similarities and differences between both approaches to solving that yield some 
statements and some open questions of interest. In this sense, we have tried to 
write down as many comments as possible to clarify (as much as we can) the 
relations between both approaches to solving. 

Moreover, we have tried to put both approaches under the same hypothesis. This 
means, that our input system of multivariate polynomials F := (/i, . . . , fs) is 
well-suited for the application of either Kronecker's or Newton's approach to 
solving. Therefore we will assume the following hypotheses : 

i) The number of equations equals the number of variables (i.e. s — n above) 

a) The variety V{fi,...,fn) is zero-dimensional and co ntains exactly D 
points, i.e. the degree of V{fi, . . . , /„) (in the sense of |EIci83|]) is exactly 
D. 

Hi) The iiT— rational points in Vk(/i, ■ ■ • , fn) are smooth with respect to the 
system F := (/i, . . . , /„), i.e. for every C G l^if (/i, • ■ • , /n), the jacobian 
matrix 

mo f 1^(0 

/ l<i,j<n 

is a non-singular matrix DF{C,) E GL{n, K). 

iv) The sequence /i, . . . , /„ e ...,/„] is a reduced regular sequence, 

i.e. for every i, 1 < i < n — 1, the ideals (/i, . . . , fi) are radical ideals of 
codimension i in Q[Xi, . . . , X„]. 



v) The degrees of the input polynomials satisfy deg(/i) < 2, for 1 < i < n. 

It must be said that constraints i) and iv) are not relevant for Kronecker's 
approach to solving . Apply i ng the iter ative v ersion of Bertini's Theorem (as 
described in [Mor97, Hag98, HMPSj or |GS9£]) we can easily reduce the over- 
determined input system to a system satistymg properties i) and iv). Anyway, 
we prefer to keep this hypothesis to simplify exposition, notations, and - hope- 
fully - reading. 

The rest of the introduction presents the new results, classified into three main 
categories : 

i) Newton's approach to solving. 

Here, we show how to extend t he approximate zero theor y introduced 



by S. Smalc in | Sma81 | (cf. also 
devclqppcd in colabqration with 
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ii) Kronecker's approach to solving. 

This recalls Kronecker's approach to solving and shows the main state- 
ments which relate both approaches by means of an algorithm based on 
the (or LLL) r eduction procedure (as introduced in |LLL82| and used 



iii) 



in |KLL84, Lcn84|) 



Application: Computation of splitting field and Lagrange resolvent. 

Finally we exhibit an algorithm that combines both approaches to com- 
pute efficiently the splitting field of an univariate polynomial equation and 
also the corresponding Lagrange resolvent. 



2 Newton's approach to solving 

Let Mir b e a proper class of absolute values on the number field K in the sense 



of | Lan83 |. For every v G Mk we have an absolute value \ ■ \i, . K 
The class A/k is chosen such that it satifies Weil's product formula with respect 
to well-defined multiplicities. We denote by 5 C Mk the set of sub-indices 
V e Mk such that the absolute value | • 1^, is archimedean and, consequently, by 
Mk \ S the class of sub-indices v G Mk such that | • |^ is non-archimedean. 
For every v S Mk, we shall denote by K^, the completion of K with respect to 
the absolute value | • We also denote by | • |i, : — > R the corresponding 
extension of | • j^^ to the completion K^. 

Let C g Vft:(/i: ■ • ■ I fn) be a smooth iiT— rational point of the zero-dimensional 
complex algebraic variety V{fi, . . . , /„). We are interested in approximating C, 
using iterations of the Newton operator. Therefore, we introduce the Newton 
operator of system F as the following list of rational mappings : 




fi{Xi, . . . , X„) ' 

> fn {Xi , . . . , Xn ) / 

An approximate zero z in A"" for the system F with associate zero C e Vr-(/i, . . . , /„) 
with respect to the absolute value | • j^^ is a point such that the sequence of it- 
erates of the Newton operator is well-defined and converges quadratically to C. 
Roughly speaking, an approximate zero z G if" with associate zero C £ K" is 
a point which lies in the basin of attraction of the actual zero ^ with respect to 
the Newton operator Np- Formally, we define approximate zeros as follows : 

Definition 1 Let F := (/i,...,/„) be a system of multivariate polynomials 
with integer coefficients : fi G Z[Jsri, . . . , X„] for 1 < i < n. Let v € Mk 

define an absolute value \-\^ . K > R. Let C, £ Vft:(/i, . . . , /„) be a smooth 

K— rational point (i.e. DF{(^) G GL{n,K)). Let z := (zi,...,z„) G be 
an affine point. We say that z is an approximate zero of the system F with 
associate zero C G with respect to the absolute value \ ■ \^^, if the following 
properties hold : 



DF{z) G GL{n, K) is a non-singular matrix. 
The following sequence is well-defined : 

zi := Np{z) G K", and Zk := iV_F(zfe-i) for k < 2. 

For every fc G N, k > 1, the following inequality holds : 

\\zk-C\\i, < -^^\\z-CL, 

where \\ ■ Wu : K^, > R is the corresponding norm associated to the 

absolute value I • I,,. 
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From a computational point of view, we want to compute approximate zeros of 
smooth if —rational points and we want to write them over a finite alphabet. 
In particular for every smooth if— rational zero G (/i, • ■ • , fn) ^^nd every 
absolute value ly S Mk, we consider a subfield L of if , such that the completion 
L^, of L with respect to the absolute value | • j^^ contains the entries of ^, namely 
C G ij". Thus, we look for approximate zeros z e with associate zero C G ij". 
Let us observe that if the absolute value \ ■ \^ is archimedean, we may fix L to 
be L := Q[i\. Moreover, we are interested in the heights of approximate zeros 
z G L"^ with actual zeros C G In the case where L = Q[i], the height of a 
point z e Q[i]" essentially equals its bit length (i.e. the number of tape cells 
in a Turing machine required to write down the number list of digits describing 
z). In the sequel, we shall therefore identify the logarithmic height ht{z) and 
its bit length. 

A first relevant task consists in stating conditions which are sufficient for veri- 
fying the property of being an approximate zero. This is achieved by means of 
a local condition based on a quantity (called 7) , which is essentially vi elded bv 



the Lipschitz constant appearing in the inverse mapping Theorem (cf. [Dem89|, 
Ch. 1, for i nstance ). These ideas were introduced by S. Smale in the early 



eighties (cf. [Sma81|) and deeply developped in the series of papers written by 
M. Shub and S. Smalc |SS8| to |SS94bt ) below). 



With the same notations as above, let v G Mk be an absolute value on the field 
if. We define the quality 7 : 



lu{FX) ■■= sup 

fc>2 



{DF{C))-\D^^)F{0) 



k\ 



where the norm is considered as the norm with respect to the absolute value 
\-\^ of the multilinear operator 

DF{C)~^D^^^F{0 : {K^ f — > if^'. 

This quantity yields a locally sufficient condition for having an approximate 
zero. This statement is known as the 7— Theorem and it holds equally true for 
archimedean and non-archimedean absolute values. 

Theorem 2 (7— Theorem) With the same notations and assumptions as be- 
fore, let F := (/i,...,/n) be a sequence of multivariate polynomials with co- 
efficients in if. Let ^ G Vr-(/i, . . . , /„) be a smooth K^rational zero (i.e. 

DF{C,) G GL{n, if) is a non-singular matrix). Let \-\^ : if > be an 

absolute value on if. For every z G if" satisfying the inequality : 

|lC-z||.7.(i^,C)<^^ 

holds : z is an approximate zero of the system F with associate zero C, with 
respect to the absolute value | • 1^. 

The proof of this statement fo llows ste p by step the proof of the usual 7— Theorems 



(cf. the compiled version in [BCSS98|) 



To establish upper and lower bounds for the bit length of approximate zeros, we 
have established several technical statements. One of them is an extensi on to t he 



non-archimedean case of the well-known Eckardt & Young Theorem [EY36| 



Let V G Mk be an absolute value over if and if^ the completion of if with 
respect to the absolute value | • |^,. Let us denote by C Mn{Kv) the variety 
of singular n x n matrices with entries in if^. Similarly, let S be the subset of 
of all singular n x n matrices with entries in if. Finally, let 

) : MniK,) X Mn{K,) R+ 

be the Frobenius (also called Hilbert-Weil) metric on A^„(if^) with respect to 
the absolute value | • | . Then, the following Theorem holds : 
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Theorem 3 (Eckardt &l Young) Let v £ Mk he an absolute value. For ev- 
ery non-singular n x n matrix A G GL(n, K), the following equality holds : 

Sf^A, S) = di^HA, E,) = inf{df M) : M e E} = 

For every multivariate polynomial / e . . . ,X„] with integer coefRcients, 

we define its logarithmic height ht{f) as the logarithm of the maximum of the 
absolute values of its coefficients. This notion introduced, we have the following 
statement which shows lower bounds for the bit length of approximate zeros. 

Theorem 4 (Lower Bounds) Let fi, ■ ■ ■ , fn G '^[Xi, . . . , X„] be a sequence 
of multivariate polynomials. Let us assume that the following properties hold : 

i) max{deg(/i) : 1 < « < = 2 , 

ii) ht{fi) < h for 1 < i < n. 

Let C e yR-(/i, . . . , /„) be a smooth K— rational point of the system F := 
(/i,...,/„). Let \ -\ . K — > be an absolute value defined on K, and 
let L CI K be a number field such that ^ G L". Then, for every z G L", z ^ C, 
satisfying : 

lk-CII.7.(i^,C)<^^ 

the following inequality holds : 

ht{z)> ^^^^ (log7.(f,C)-[£ : m (5 log n + 2/1) -3). 

Using Theorem above, the following inequality also holds : 

ht{z) > ^^^\ (log4^)(Jf (C)-\I].) - [L : Q](71ogn + 3/i)-5). 

Moreover, in the case where L — Q[i] is the field of Gaussian rationals, the two 
previous lower bounds may be rewritten as : 

ht(z) > - (logj^(F,C) ~ (101ogn + 4/1 + 3)), and 
6 

ht{z) > i (logdf)(i?F(C)-\E.) - (141ogn + 6/i + 5)) . 

Let us observe that the "negative terms" in the previous lower bounds are linear 
in the input length (i.e. the bit length of the input system F := (/i, . . . , /„)) 
whereas the "positive part" depends semantically on the smooth ii'— rational 
solution C G . . . , /„). 

Last, but not least, we may also show a few lower bounds for the average height 
of approximate zeros associated to a Q— definable irreducible component of the 
solution variety V{fi, . . . , /„). To this end, we introduce some additional nota- 
tions. Let /i, . . . , /„ G . . . , Xn] be a sequence of multivariate polynomials 
satisfying the hypotheses i) to v) above. Let ( G K" be a smooth if— rational 
zero of the system F := (/i, . . . , /„). Let V V{fi, ...,/„) C C" be the 
algebraic variety given as the common zeros of the polynomials /i, . . . , /„. Let 
V(; C 1/ be the Q— definable irreducible component of V that contains C • Let 
us assume D := degiVc^) be the number of points in V(j. Let us observe that 
D = deg{() < [K : Q]. Let us assume 

Vc :={Ci,...,Cd}. 
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Let II • II : — > M be the standard hermitian norm induced in /C" by the 
inclusion i : K ^ C A sequence of points z := (zi, . . . , zd) G Q[«]"^ is said to 
be an approximate zero of the system F with associate variety that satisfies 
the 7— Theorem if for every i, 1 < i < D, the foUowing holds : 

II rll <- 



where 7(i^, Ci) is the quantity associate to the hermitian norm || • ||. 
For every given approximate zero z :— (zi, . . . , zd) <= of the system F 

with associate variety V"^, the average height (also the average bit length) of z 
is defined in the following terms 

1 ^ 

htaviz) X! '^*(^»)- 

1=1 

Finally, let us denote by Zk C K the ring of algebraic integers of the number 
field K. Then, we have the following lower bound for the average bit length of 
approximate zeros with associate variety : 

Proposition 5 With the previous notations, let C. G VK{fi, ■ ■ • , fn) be a smooth 
K— rational with entries in Ijk, *-e- C S Z^. Let us also assume that for every 
archimedean absolute value \-\^ (i. e. v £ S), the following holds : 

3||C|U7.(i^,C) >3-%/7. 

Then the average height of any approximate zero z € Ql?]"^ of the system F 
with associate variety satisfies the following inequality : 



htav{z) > ^ 



M(C)-(ilogn + log2) 



In order to illustrate the meaning of this lower bound, we give here a few Corol- 
laries. 

Corollary 6 With the same notations as in Proposition aftoue, let ( £ TLk H 
^Kifi, ■ ■ ■ , fn) be a smooth K— rational zero of the system F : 



let us assume that for every archimedean absolute value 
for every v £ S), the following holds : 



ifl, 

: K 



,/„) and 
(i. e. 



Then the average height of any approximate zero z £ Q[«]" of the system F 
with associate variety V"^ satisfies the following inequality : 



htav{z) > 



htio 



(-logn + log2) 



Moreover, the previous techniques show how to deform a given system of mul- 
tivariate polynomials by means of a single additional equation of low degree in 
such a way that the average bit length of the new system is essentially greater 
than the height of the actual zero you want to approximate. 

Corollary 7 Let F (/i,...,/„) be a system of multivariate polynomials 
with integer coefficients satisfying the conditions i) to v) given on page [J. Let 
C = (Ci, . . . , C»i) G '^K n Vff (/i, . . . , fn) be a smooth K— rational zero whose 
coordinates are algebraic integers. Let us now define the system of polynomial 
equations in n + \ variables : 



G:= (5i,...,5„+i) e (Z[Xi, 
given by the following rules : 



n+1 
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• 9i ■= ft e "^[Xi, . . . , Xn+i] for every i, I <i <n, 

• gn+l '■= {Xn+1 — Xn) (X„+i — (X„ + 1)) . 

Let C G VKigi, ■ ■ ■ , gn+l) n be the affine point given by : 

C (Cli ■ • ■ J CniCn) G 

Let Vq' C V{gi, . . . , (7„+i) &e t/ie Q—definable irreducible component ofV{gi, . . . , gn+i) 
containing C'. Then, the average height of any approximate zero z G Q[z]("+"'^^^ 
of the system F with associate variety Vqi satisfies the following inequality : 



htav{z) > ^ 



ht{C)-{^\og{n + l) + \og2) 



We can construct several examples where all of the previous lower bounds for 
the bit length of approximate zeros apply. For instance we give here an exa mple 
inspired by a classical univariate example due to M. Mignotte (cf. [Mig89|) : 

Example 1 (Using log 7 as in Theorem |^) Let us consider the system of 
multivariate polynomials F := (/i, . . . , fn+i) given by the following rules : 

. /i := Xi - 2, 

• fi :~ Xi — Xf_i for every 2 < i < n — 1, 

• fn ■= Xn-^-l — X^, 



• /n+1 Xn+iXn — 2(X„_iX,i — 1 



|2 



This system F has three solutions in C"^"'^, where two of them, say Ci, C2 G M"^"'^, 
satisfy the following inequality : 

||Cl-C2||<-|r^<rJrT- 



Thus, using the well-known relation between the quantity 7 and the separation 
of roots, we may conclude : 

2 >„ 3-V7 

j^rrr > IIC1-C2II > 



By Theorem^, we conclude that for all approximate zeros Zi,Z2 G Q[i]"^"'^ of 
the system F^ associated to ^1 , ^2 respectively and satisfying the corresponding 
7— Theorem, the following holds : 

ht{zi) > i (2"-i - 21og(n + 1)) - 0(1). 

Therefore this example (an several others wich should be constructed) allows us 
to conclude the following Corollaries ^ to |l0| . 

Corollary 8 Computing approximate zeros in Q[i] for archimedean absolute 
values, using binary encoding of the output requires exponential running time 
and exponential output length, and these two lower bounds cannot be improved 
with this encoding. Namely, computing approximate zeros with binary encoding 
is in the complexity class EXTIME \ P. 



Corollary 9 Floating point encoding of approximate zeros requires exponential 
number of digits and this lower bound cannot be improved. Namely, floating 
point encoding does not suffice to compute approximate zeros of systems of mul- 
tivariate polynomial equations. 
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As continuous fraction encoding of numbers in Q[i\ is close to the binary encod- 
ing, we easily conclude the following : 

Corollary 10 Computing approximate zeros in Q[i] for archimedean absolute 
values, using continuous fraction encoding of the output requires exponential 
running time and exponential output length, and these two lower bounds cannot 
be improved with this encoding. Namely, computing approximate zeros with 
continuous fraction encoding is in the complexity class EXTIME \ P. 

These lower bounds suggest that a central point of interest should be to study 
the bit length of approximate zeros satisfying the 7— Theorem. In order to shed 
some light in this direction, we prove the following statements : 

Theorem 11 (Upper Bounds) Let fi, ■ ■ ■ , fn G '^[Xi, . . . , Xn] be polynomi- 
als with integer coefficients. Let us assume that the following properties hold : 

• max{deg(/i) '. 1 <i < n} <2, and 

• ht{fi) < h for 1 < i < n. 

Let C £ Vr-(/i, ...,/„) be a smooth K— rational point. Let | • |^ : K — > R-|_ be 
an absolute value on K . Then, the following inequality holds : 

log7^(F,C) < 3[is: : Q]n(n2 + 41og7i + ;i + /it(C)+3). 

In particular, we show the following estimate for the bit length of approximate 
zeros in (Q[i]" : 



Corollary 12 (Upper bound on the bit length of approximate zeros) 

With the same assumptions and notations as in Theorem ^ above, let ^ G 
Vxifi, ■ ■ ■ 1 fn) be a smooth K—rational zero, and let | • |^ be an absolute value 
on K. Let L i- K be a number field such that € L". Then there exist approxi- 
mate zeros z e L" of the system F :— (/i, . . . , /„) with approximate zero C with 
respect to the absolute value \ ■ \^, such that the logarithmic height ht{z) of z is 
at most linear in the following quantities : 



1 



[L : 



■log|AL| + [X 



(n^ + h + nhtiO) , 



where jA^j is the absolute value of the discriminant of the field L. 

Moreover, in the case where L = Q[i] (for instance, if \ ■ \^ is archimedean), 
there exist approximate zeros z € Q[i]" for the system F with associate zero ^ 
with respect to | • j^^ such that their bit length is at most linear in the following 
quantity : 

[K : Q]n (n^ + h + nht{()) , in other words : 
ht{z) <0{[K : Q]n {n^ + h + nht{C))) . 
Let us observe, that these two upper bounds above (i.e. Theorem ^ and Corol- 



lary 
whic 



I2D depend mainly on the input length n h 



and on two parameters 



1 m turn depend on the actual zero to approximate : the degree [K : Q] of a 
field containing the coordinates of the zero and the logarithmic height of the zero 
ht{C,). These t wo qua ntiti es arc bounde d respectively by the geometric Bezout 
inequ alitv fcf . |Hci83| or |Ful84 . Vog84 n and the arithmetic Bezout inequality 
(cf. pGS93| , pluyil , ]Phigl^ ' pEi9^ or pCP9^ ^P9|, |Som9^ , lHag98| , |imPS[ , 
for instancej. Moreover, conibinmg these two upper bounds ( Iheorern and 
Corollary 12) with the previously shown lower bounds and several examples 
may concluc 
are optimal. 



3) with the previously shown lower bounds and several examples, we 
that the upper bounds shown in Theorem |ll] and Corollary ^ 



On the other hand, the 7— Theorem above has some aesthetic consequences 
which we may explain in terms of the existence of a universal radius of conver- 
gence independent of the aboslute value under consideration. To this end, we 
recall the well-known Implicit Function Theorem for complete noetherian local 
rings in the following terms : 
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Theorem 13 (Non archimedean Basin of Attraction) Let F :— (/i, .... /„) 
Z,[Xi, . . . , Xr,]" be a system of multivariate polynomials satisfies the hypotheses 
of TheoremiTi above. Let v S Mj^ define a non-archimedan absolute value \ ■ \^ 
on K . Let us also assume that the restriction 



defines a p~adic absolute value, where p € N is a prime number. Let C G K"' 
be a smooth K— rational zero of the system which lies in the closed unit shere 
of K'\ I.e. 

CeB,(0,l):={a;G/^" : < 1}. 

Let us finally assume that \ dei DF{C,)\i, = 1. Then, for every z e i?i,(0, 1) 
satisfying 

p 

holds : z is an approximate zero of the system F with associate zero Q with 
respect to the absolute value | ■ 

This statemen t is nothing but the usual Hensel Lemma in local algebra (of. 
[|ZS5S , Mor97 |, for instance). However, this statement has a drawback : The 
radius ot the basin of attraction centered at C depends on the concrete absolute 
value I • Ij,. The 7-Theorem [TT| above shows that there exists a universal radius, 
which depends only on the system F and the smooth if— rational zero, but does 
not depend on the absolute value. 



In order to prove this claim, let us introduce the following quantity ^{F,(^). 
With the same notations and assumptions as above, let us define the universal 
quantity 

C\ WW] 
Y[ max{l,7.(i^,C)r-' 

Let lis observe, that this quantity is well-defined and finite according to Theo- 
rem |ll| above. Moreover, it does not depend on any particular absolute value 
under consideration. Thus, we may conclude the following Theorem : 

Corollary 14 (Universal 7— Theorem) With the same notations and assump- 
tions as in Theorem^ for every z G Q[i]" and every absolute value \ ■ \^ satis- 
fying the following inequality 

\\z-C\\.liF,0<^^ 

holds : z is an approximate zero for the system F with associate zero C and with 
respect to the absolute value v € Mk. 

Let us point out that the existence of such a universal quantity does not imply 
the existence of a universal basin of attraction independent of the absolute value 
under consideration. In fact, we show the following (expected) statement : 

Corollary 15 Let F := (/i, •..,/«) be a sequence of multivariate polynomi- 
als with integer coefficients satisfying conditions i) to v) of page^. Let C G 
^A'(/i7 • ■ • , fn) be a smooth K— rational zero. The only point z G that sat- 
isfies the universal j— Theorem near 7 for all absolute values in Mk is z = Q. 
Namely, for every z G if" satisfying the following inequality for every v £ Mk 

II ^ 3-x/7 

holds z ~ 
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3 Kronecker's approach to solving 



In |Kro82|, Kronecker introduced a notion of solution of unmixed complex 
algebraic varieties, which we are going to reproduce here. Let fi, ■ ■ ■ , fi G 
Z[Xi, . . . , Xn] be a sequence of polynomials defining a radical ideal (/i, . . . , fi) 
of codimension i in C[Xi, . . . , Xn]- Let V := V{fi, . . . , fi) C C" be the com- 
plex algebraic variety of codimension i given by the common zeros of the fi. 
A solution of is a birational isomorphism of V with some complex algebraic 
hypersurface in a space of adecuate dimension. 

Technically, this is is expressed as follows. First of all, let us assume that the 
variables Xi, . . . , Xn are in Noether position with respect to the variety V, i.e. 
we assume that the following is an integral ring extension : 

Q[Xi, . . .,Xn-^] Q[Xi, . . . , Xn]/ifu ■ ■ ■ : fi)- 

Let u := \n-i+iXn^+i + • • • + A„X„ G Q[^i, • ■ • , Xn] be a linear form in the 
dependent variables {Xn^i+l^ ■ ■ ■ i ^n}- Thus we have a linear projection 

: C" — > C""'+i 

(xi , . . . , Xn) ' ^ (^1 , ■ ■ ■ , Xn—i , , . . . , Xn)) • 

Let us also consider the restriction 

U\v -V - 



The linear form u is called a primitive element if and only if the projection hi \v 
defined a birational isomorphism of V with some complex hypersurface in 
i£n-i+i -^j^j^ minimal equation Xu S 'Q.l^i, ■ ■ ■ TXn-i,T]. Then, a Kronecker 
solution of variety is a description of the primitive element m, the hypersurface 
Hu through the minimal equation Xu^ and a description of the inverse of the 
birational isomorphism, i.e. {lA Formally, this list of items may be 

described as follows : 

• The list of variables in Noether position which includes a 
description of the dimension of V . 

• The primitive element u := Xn-i+iXm+i + • • • + A„X„ given by its coef- 
ficients in Z. 

• The minimal equation of the hypersurface namely 

Xu G '^[Xi, . . . , Xn-i, T]. 

• A description of (U \v)^^- This description is given by the following list : 

— A non-zero polynomial p € Z[Xi, . . . , Xn~i]. 

— A list of polynomials Vj € Z[Xi, . . . , Xn-i, T], n — i + I < j < n, 
such that the degrees with respect to variable T satisfy degrp{vj) < 
degT(Xu), for every j,n-i + l<j<n. 

such that the following holds 

{U \v)^^ix,t) := (xi, . . . ,Xn-i,P^^ix)Vn-i+l{x,t), . . . , p^'^ {x)Vnix , t)) , 

where x {xi, . . . ,Xn-i) G C"^* and t G C. 

Kronecker conceived an iterative procedure to solve multivariate systems of 
equations F := (/i, . . . , /„) defining zero-dimensional complex varieties, which 
can be described in the following terms : 

First, you start with system (/i) and you "solve" the unmixed variety of codi- 
mension 1, V{fi) C C". Then you proceed iteratively : From Kronecker's 
solution of the variety V{fi, . . . , fi) you eliminate the new equation fi+i to ob- 
tain a Kronecker solution of the "next" variety V{fi, . . . , fi+i)- Proceed until 
you reach i = n. This iterative procedure has two main drawbacks which can 
be explained in the following terms : 
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First of all, the space problem arising with the representation of the in- 
termediate polynomials. The polynomials %„, p and Vj are polynomials 

of high degree (eventually of degree 2*) involving several variables. Thus, 
to represent them, one has to handle their coefficients, which amounts to 
the following quantities 

+ n — i + 1 
n — i + 1 

which for i :— n/2 amounts to more than 2" /"^ coefficients of great bit 
length. 

• Secondly, Kronecker's iterative procedure introduced a nesting of interpo- 
lation procedures required for the iterative process and the linear change 
of cooerdinates required by every computation of the Noether normali- 
sation. This nesting of interpolation procedures is difficult to avoid and 
increases the running time complexity. 

Therefore, the procedure was forgotten by contemporary mathematicians and 
hardly mentioned in t he litera ture of algebraic geom etry. M acaulay quotes Kro- 
necker's procedure in Macl6| and so does Konig in iKon03| . But both of them 



thought, that this procedure requires too much running time to be efficient, and 
it was progressively forgotten. Traces of this procedure can be found spread 
over the algebraic geometry literature without giving to it the required rele- 
vance. For example, Kronecker's notion of solution was used by O. Zariski in 
Zar95| to define dimension of algebraic varieties, claiming that it was also used 



solving wi thout 



prcvi- 



in the same form by Severi and others. 
In 1995, two works rediscovered Kronecker's approach to 

ous knowledge of it's existing ancestors. These two works 

were able to overcome the first drawback (space problem of representation) of the 
previous methods. The technical trick was the use of a data structure coming 
from semi-numerical modelling : straight-line programs. This idea of repre- 
senting polynomials by programs that evaluate them go es b ack to former works 



GHMP95t| and Par95 | 



of the same research group (as [|GH93| , |FGS95i [KP9l or |KP96| ). Moreover 



these id eas were the natural continuation of the ideas previously developped in 



GH91 |. 



To overcome the second drawback (Nesting), the authors introduced a method 
based on Newton's method applied in a non-archimedean context (the approx- 
imate zeros in the corr esponding non-archimedean basin of attraction were 



called Lifting Fibers in [GHH+97 |). This result was obtained in the two pa- 
pers iGHH+97t |GHM+98|] . The key trick to avoid the nesting of interpolation 
procedures was based on Hensel's Lemma (also Implicit Mapping Theorem). 
Perhaps, the following statement could explain the relations existing between 
Hensel's Lemma and Approximate Zero Theory. To this end, let us introduce 
a few more notations. Let /i, . . . , € C[Xi, . . . , X„] be a sequece of polyno- 
mials defining a radical ideal of codimension r in C[Xi, . . . , X„]. Let us assume 
that the variables Xi , . . . , X„ are in Noether position with respect to the ideal 
/ :— (/i, . . . , /r), i.e. assume that the following ring extension is integral 

C[Xi,...,X„_,] -^C[Xi,...,X„]//. 

Let P := {pi, . . . ,Pn~r) G C"^'' be an affinc point, let Op be the ring of formal 
power series at P, and let A4p be the field of fractions of Op. Then, the 
following is finite ring extension 

Mp^B:= Mp[X.n-r+l, . . . , ^„]/(/l, . . . , fr), 

and i? is a zero-dimensional A^p— algebra. Thus, it has some sense to look for 
approximate zeros of the solutions in A^p of the system of polynomial equations 
F := {fi, . . . , fr). The following statement about Hensel's Lemma explains the 
connections existing between Kronecker's solving and Approximate Zero Theory. 
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Theorem 16 (Hensel's Lemma) With the same assumptions and notations 



as above, let ( Aip be a solution of the system F . 
be usual non-archimedian norm in M^p. Let <C{Xi, . 
rational functions. Then, for every z € <C{Xi, . . . , X„_ 



Let II -11 -.M] 

. ,Xn-r) be the field of 
~Y , if \\z\\ < 1, and 



lk-CII< 



1 

2' 



then z is an approximate zero for the system F :— (/i, . 
zero C e M^p. 



, fr) with associate 



Unfortunately, those two works GHH+97 , GHM+98 | introduced (for the Lifting 
Fibers) run time requirem ents wh i ch depend on the heights of th e intermediate 
varieties (in the sense of 



This draw- 

gger numbers were 

represented by straight-Une programs and the following result established : 



Phi9i 



back was finally overcome m the paper | GHMP97 1 



[Phi94 ^ifi9^ ^ ^om98 
I1MP97| where mtegei 



Theorem 17 ; GHMP9^ There exists a bounded error probability Turing ma- 
chine M which performs the following task : Given a system of multivariate 
polynomial equations F := (/i, . . . , fn), satisfying the following properties 

• deg(/i) < 2 and ht{fi) < h for I < i < n, 



• the ideals (/i, . 
for 1 < i < n - 



1, 



, fi) are radical ideals of codimension i in the ring 



• the variety V{fi, . . . , fn) C C" is a zero-dimensional complex algebraic 
variety. 

Then, the machine M outputs a Kronecker solution of the variety V{ fi, . . . , /„). 
The running time of the machine M is polynomial in the following quantities 

5{F)nh, 

wh e.re () is the maximum of the degrees of the intermediate varieties ( in the sense 
of i HeiSjjl ), namely 



6{F) -.^ max{deg{V{fi, ...,/,)) : 1 < z < n ~ 1}. 

It must be said that the coefficients of the polynomials involved in a Kronecker 
solution of the variety V(fi, . . . , /„) are given by straight -line programs that 
evaluate integer numbers. However, the complexity estimates for the Turing 
machine M are independent from the height. 

Our attempt in these pages is to compare this approach to solving developped 
by Kronecker to that of Newton as described in the previous Subsection. 

The exposition of new results starts with a sr nall impr ovement of the Witness 
Theorem in jHSSCH and [ |BCSS96i (cf. also iBCSS98| ). When deahng with 
straight -line program data structures, some relevant technical methods of com- 
parison must be developped. These methods are known as probabilistic zero 
tests for polynomials giv en bv s traigh t -line p rograms. Examples of these tests 



are those introd uced in [3ch79 
duced in |HS8C] for the case o: 



Zip79| , |HS80[ and the Witness Theorem, intro- 
poiynomials with integer coefficients, and in 



BCSS96| lor polynomials with coefficients in a number field. 



As we already had to introduce a few technical notions and methods spread over 
the literature of number theory, numerical analysis, algebraic complexity theory 
and elimination theory, we can give without introducing additional material for 
free the following improvement of the estimates for the Witness Theorem. 

Theorem 18 (Witness Theorem) Let P e K[Xi, . . . , Xr^ be a non-zero 
polynomial evaluable by a non-scalar straight-line program F of size L, non- 
scalar depth £ and parameters in J-' Q K . Let luq Cz K be such that the following 
holds : 

ht{uJo) > niax{log2,/it(J^)}. 
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Let N £ N be a non-negative integer such that 

log, N > \og,{e +!) + {£ + 2) (log2 log2(4L)) . 

Let us define recursively the following sequence of algebraic numbers (known as 
Kronecker's scheme) : 

uji :— lUq , and for 2 < i < n, uJi := ujf^i- 

Then, the point lo :— (wi, . . . ,ujn) G is a witness for P (i.e. P{uJ_) ^ 0). 

Moreover, we observe that approximate zeros are succint encodings of generic 
points of the variety V{fi ,...,/„). This means that for every smooth ii'— rational 
zero C G (/i, . . . , /„), the binary encoding of an approximate zero z £ Q[i]" is 
sufficient information to compute the Q— irreducible component of V{fi, . . . , /„) 
containing (.. In more precise terms we show the following statement : 

Theorem 19 (From Approximate Zeros to Geometric Solution) With the 
same assumptions as in Theorem [7J' above, there exists a bounded error prob- 
ability Turing machine M , such that taking as input the binary encoding of an 
approximate zero z G Q[i] of the system F with associate zero Q £ Vxifi, . . . , fn) 
for an archimedean absolute value \-\^ (where v £ S), M outputs a Kro- 
necker solution of the Q— irreducible component W of V{fi, ...,/„) containing 
Moreover, the running time of this probabilistic Turing machine is polynomial 
in the following quantities 

deg{W) {n h ht{z)ht{0) , 
where deg(W^) is the degree of the irreducible component W containing C,. 

The key idea for the proof of this Theorem is the use of the L^ (or LLL) 
reduction algorithm. 

Conversely, as approximate zeros may depend on the height of the actual zero 
they approximate, we could be interested in the computation of approximate 
zeros for actual zeros with small (bounded) height. This is done in the following 
statement. 

Theorem 20 (Prom Kronecker's solution to Newton's solution) There ex- 
ists a bounded error probability Turing mahcine M which performs the following 
task : Given a sequence of polynomial equations F := (/i, . . . ,/„) of degree at 
most 2 and height at most h, and given a positive integer number H £ N in 
binary encoding, the machine M outputs approximate zeros for the archimedean 
absolute value \ ■ \ : K — > M induced on K by the inclusion i : K ^ C for all 
those zeros C £ Vxifi, . • . , fn), whose logarithmic height is at most H , i.e. 

htiO < H. 

The running time of M is polynomial in the following quantities : 

{ndd{F)) + {Dnh H), 

where the notations are the same as in Theorem ^ before. 

A proof of this statement is based again on an application of the L^ reduction 
algorithm. 
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4 Application : Computation of splitting field 
and Lagrange resolvent 

Combining both Kronecker's and Newton's approach to solving, we exhibit an 
efficient procedure for computing the splitting field and the Lagrange resolvent 
of an irreducible monic univariate polynomial / G Q[X] of degree d. Let us 
recall that the splitting field of / is the minimal number field K{f) containing 
the field of rational numbers Q and all complex roots of / (i.e. the minimal 
number field where / splits completely, also called the normal closiire of the 
equation / = 0). This normal closure /!'(/) is nothing but the Galois field 
of / and it satisfies where GalQ{f) is the Galois group of the polynomial /. 
The splitting field of / can be indentified with an irreducible component of the 
zero-dimensional algebra (known as the universal decomposition algebra) 

A := Q[Xi, Xd\/{ao - ao, ■ ■ ■ , <Jd-i - Ud-i), 

where ctq, . . . , Cd-i are the elementary symmetric functions and / is written 

as f{X) := oo + aiX + • • • + a^^iX"^^^ + X'^. Let us also observe that the 
Lagrange resolvent is nothing but the Chow (or Cayley) elimination polynomial 
of the zero dimensional resiude algebra A/m, where m is a well-chosen maximal 
ideal of A. Therefore, we can also show the following Theorem as a consequence 
of the comparison between Newton's and Kronecker's approach to solving : 

Theorem 21 (Splitting Field and Lagrange Resolvent) There exists a prob- 
abilistic Turing machine, which for every given univariate polynomial f G 
of degree at most d and logarithmic height at most h computes the following 
items : 

i) Approximate zeros in Q[i] of all zeros of f, 

a) a geometric description of the splitting field K{f) of the polynomial f , 
Hi) and the Lagrange resolvent of the equation / = 0. 
The running time of M is polynomial in the following quantities : 

UGalif)){dh). 
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